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ABSTRACT 

We quantify the accuracy with which the cosmological parameters characterizing the 
energy density of matter (f2 m ), the amplitude of the power spectrum of matter fluctu- 
ations (erg), the energy density of neutrinos (f2„) and the dark energy equation of state 
(wo) can be constrained using data from large galaxy redshift surveys. We advocate a 
joint analysis of the abundance of galaxies, galaxy clustering, and the galaxy-galaxy 
weak lensing signal in order to simultaneously constrain the halo occupation statistics 
(i.e., galaxy bias) and the cosmological parameters of interest. We parameterize the 
halo occupation distribution of galaxies in terms of the conditional luminosity func- 
tion and use the analytical framework of the halo model described in our companion 
paper (van den Bosch et al. 2012), to predict the relevant observablcs. By perform- 
ing a Fisher matrix analysis, we show that a joint analysis of these observables, even 
with the precision with which they are currently measured from the Sloan Digital Sky 
Survey, can be used to obtain tight constraints on the cosmological parameters, fully 
marginalized over uncertainties in galaxy bias. We demonstrate that the cosmological 
constraints from such an analysis are nearly uncorrclated with the halo occupation 
distribution constraints, thus, minimizing the systematic impact of any imperfections 
in modeling the halo occupation statistics on the cosmological constraints. In fact, 
we demonstrate that the constraints from such an analysis are both complementary 
to and competitive with existing constraints on these parameters from a number of 
other techniques, such as cluster abundances, cosmic shear and/or baryon acoustic 
oscillations, thus paving the way to test the concordance cosmological model. 

Key words: galaxies: cosmology — galaxies: halos — galaxies: structure — dark 
matter — methods: statistical 



1 INTRODUCTION 

Of the various cosmological models in the literature, the 
ACDM model has withstood a large variety of observational 
tests which have grown more and more stringent over time. 
According to this model, at present times, dark energy and 
dark matter form the dominant component of the energy 
density budget of the Universe. Ordinary matter forms a 
sub-dominant component and is primarily present in the 
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form of gas in the intergalactic medium and around galaxies. 
Obtaining precise constraints on the energy densities of the 
various components of the Universe (quantified by the en- 
ergy density parameters for matter [fl m ], dark energy [Qa], 
baryonic matter [fib] and neutrinos [fii/]), is of great im- 
portance, in order to understand the expansion history and 
future fate of the Universe. 

Any successful cosmological model also requires to ex- 
plain the formation of structure in the Universe. The ACDM 
model assumes the presence of nearly scale invariant, tiny 
initial fluctuations in the matter density, presumably gener- 
ated during inflation. Structure grows hierarchically in this 
model with small scales collapsing first followed by larger 
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and larger scales as the Universe grows old. The precise de- 
scription of the inhomogeneous Universe requires the knowl- 
edge of the index of the power spectrum of initial fluctua- 
tions, 7i s and its amplitude (quantified in terms of the pa- 
rameter as). 

Observations of the fluctuations in the cosmic mi- 
crowave background on both large (e.g.. ISpergel et al] 20071; 



Komatsu et al.l 2011) and small scales (e.g., iLueker et al.l 
20101 ; iDunklev et ai]l201ll ) , of the dimmin g of distant super- 



novae as a function of t heir redshift (e.g..lRiess et al.l 19981; 
Perlmutter et al. I ll999l ; iKowalski etai] |200S| ; iKessler et al.l 



20091 ; iGuv et al.l l2010l). of the large scale structure tracec 



by galaxies (e.g.. lTegmark et al.ll20o3 ; ISwanson et alj|2010l) . 
of the abundance of massive clusters and its evolution (e.g. , 
Vikhlinin et~aH |2009| ; iMantz et all l20ld ; iRozo et all l20ld ; 



Benson et al.ll201ll ; ISehgaT et al.l l201lh . of the angular ex 



tent of the bary onic acoustic oscillation scale as a function 
of redshift (e.g.. lEisenstein et al.ll2005l; [Percival et al.ll2007l . 



|2010| ; iBlake et alj|201ll ; lAnderson et al.ll2012l ), and of the 

weak distortion of galaxies due to intervening structure (e.g 



Massev et al.1 1200 



galaxies due to intervening structure (e.g., 
3; ISchrabback et all |2010| ; iLin et~ai] 1 20 lit 



Huff et al.ll201ll ) have all contributed to enable precise con 



straints on the parameters that describe the ACDM model. 
It is remarkable that a single model with a small number 
of parameters is able to self consistently explain all these 
observations. 

Large spectroscopic galaxy reds hift surveys such a s the 
2-degree field galaxy redshift survey (|Colless et al.l200lh and 
the S l oan digital sky s urvey (hereafter SDSS, York et all 
|2000| ; lAbazaiian et al.1 120091 ) have played an important 
role of mapping out the three-dimensional structure of 
galaxies in the Universe in exquisite detail. This has al- 
lowed precise measurements of the abundance of galax- 
ies and their clustering on non-linear sc ales as a function 
of galaxy propertie s such as luminosity |Wang et al. 20081; 



01 galaxy properties such as luminosity ( W 
Zehavi et ail l201ll ; iBlanton et ail l2003bT 



Z.enavi et al.l I^Uill ; imanton et aU I2UU3P); iNorberg et al.l 
200 if), stellar mass dLi et al.1 120071; iBaldrv et al.1 1 20081 



Cole et alj200lh. or co lo ur dNorberg et alj|2002 



20081 ; ISwanson et al.1 120081 ; IZehavi et al 



Wang et al.l 



2011 



The 

halo model has traditionally been used to obtain the 
halo occupation distribution of galaxies given a cos- 
mological model based on these observations (e.g. 



galaxy clustering data by adjusting the halo occupation dis- 
tribution of galaxies. However, we also showed that these 
models make vastly different predictions for the amount 
of baryons that should occupy dark matter haloes, quanti- 
fied by the mass-to-light ratio (see also van den Bosch et al.l 
12001 120071 ; iTinker et al] 120051 ; lYoo et al.1 l200rj |. Galaxy^ 
galaxy lensing, which measures the projected galaxy-matter 
correlation function is an excellent probe of the mass-to- 
light ratios. Therefore, it was suggested that a joint analysis 
of the abundance, clustering and lensing of galaxies can be 
used to obtain constraints on cosmological parameters such 
as f2 m and erg. 

This paper is the second in a series of three papers in 
which we develop this idea further. In van den Bosch et al. 
(2012, hereafter Paper I) we describe our model for the halo 
occupation distribution of galaxies and present an analytical 
framework for the computation of the abundance of galax- 
ies, their clustering and the galaxy-matter cross correlation. 
We make extensive use of mock catalogs to validate our an- 
alytical prescription and show that our analytical model is 
accurate to better than 10 (in most cases 5) percent in re- 
producing the 3-dimensional galaxy-galaxy correlation and 
galaxy-matter correlation. These correlations can then be 
projected to predict the galaxy clustering and the galaxy- 
galaxy lensing signal measured from SDSS data. In this pa- 
per, we perform a Fisher matrix analysis to forecast the ac- 
curacy with which various cosmological parameters can be 
constrained with current data. The analysis allows us to un- 
derstand the strength of each of our data-sets as well as the 
various degeneracies between our model parameters. In Cac- 
ciato et al. (2012, hereafter Paper III), we apply our method 
to data from the SDSS, simultaneously constraining galaxy 
bias and cosmologjQ. 

This paper is organized as follows. In section [2] we 
briefly describe the data and our analytical model. In section 
[3] we describe the Fisher information matrix corresponding 
to the observables. In section [4] we present the constraints 
on our model parameters that are statistically achievable 
given the accuracy of the data. Finally, we summarize our 
results in Section [5] and provide a future outlook. 



Jing et al.ll 19981; iBerlind fc Weinberj2002l;lYang et al.ll2003l; 2 DATA AND THE ANALYTIC MODEL 



Zheng et al.l 2005; 


van den Bosch et al.l 2007 


; Zheng et al. 


2007; 


Tinker et al. 


20051; IZehavi et al. 201 ll; 


Couoon et al. 


2011; 


Leauthaud et al.ll2012h. A few studies however have 



also turned the question around to judge if these observa- 
tions can be used to constrain the cosmological parameters 



them s elves (Ivan den Bosch et al] 2003 



20051 ; lYoo et al.l 120061 ; ICacciato et al 
20121 ) 



Tinker et al 



Tinker et al 



Galaxies reside in dark matter haloes. The cosmological 
parameters determine the abundance and clustering of dark 
matter haloes of a given mass. Therefore, the observations of 
the abundance and clustering of galaxies is sensitive to var- 
ious cosmological parameters and can be used to constrain 
these parameters. However, this requires the knowledge of 
an accurate mapping between galaxies and their dark mat- 
ter haloes. In Cacciato et al. (2009), we highlighted this 
problem, by demonstrating that two different cosmologi- 
cal models, (differing primarily in f2 m and as) were able 
to simultaneously fit the galaxy abundance and the galaxy- 



Our analysis is based on galaxy observations carried out as 
part of the spectroscopic component of the SDSS. In partic- 
ular, we focus on th e accurate measurem ents of the galaxy 
luminosity function (Bl anton et al]|2003a[). the clustering o f 
galaxies in six different luminosity bins ( Zehavi et al.l 120111 ) 
on scales ranging from 0.17 /i -1 Mpc to 42.3 /i _1 Mpc, and 
the galaxy-galax y lensing signal around g alaxies in different 
luminosity bins ( Mandelbaum et al.ll2006r ) on scales ranging 
from 0.045 /i _1 Mpc to 1.81 h~ 1 Mpc. Our aim is to perform 
a Fisher matrix analysis to understand the extent to which 
these datasets constrain the cosmological parameters that 
we are interested in and the various degeneracies that enter 
our analysis. For this purpose, we choose a set of fiducial pa- 
rameters that describe the cosmology (see Table [TJ and the 



1 A preliminary version of the main results from this paper 
and Paper III w ere published in a conference proceedings by 
iMore et al.l | |2012h . 
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Table 1. Cosmological model parameters 



Table 2. CLF and nuisance model parameters 



Cosmological 
parameters 


Model A 


Model B 


Model C 


dm 


0.266 


0.266 


0.266 


n K 


o.ot 


o.ot 


o.ot 




0.02258* 


0.02258* 


0.02258* 


Civ 


0.0 1 " 


0.004 


o.ot 




0.801 


0.801 


0.801 


n a 


0.963* 


0.963* 


0.963* 


It 


0.71* 


0.71* 


0.71* 


Wo 


-l.Ot 


-l.Ot 


-1.0 



Cosmological model parameters used to generate artificial 
datasets. We use a prior given by the covariance matrix of the 
WMAP chain for the parameters marked with an asterisk. The 
parameters marked with a dagger are kept fixed when analysing 
the model for the fiducial case. We discuss variations on some of 
the assumed priors in the text. 

halo occupation distribution of the galaxies (see Table [2]). 
We use the analytical model (see Sections I2,U2.4[) to pre- 
dict galaxy abundances, galaxy clustering (in projection), 
and the galaxy-galaxy lensing signal (always using the same 
luminosity bins and radii as in the SDSS data). We assume 
that these observables are measured with a percentage ac- 
curacy that is listed in Table [3] and chosen such that the 
measurement quality is fairly similar to that of actual SDSS 
measurements. 

We perform three different analyses, that correspond to 
the ACDM model and its variants. For the first analysis, we 
assume a flat standard ACDM cosmology with a negligible 
amount of energy density in neutrinos. We primarily focus 
on the matter density parameter, fi m and the parameter 
that characterizes the amplitude of matter fluctuations, erg. 
We use priors on the secondary cosmological parameters, the 
baryon energy density parameter (fit,), the Hubble parame- 
ter (h = fzo/100 kms -1 Mpc -1 ) and the power law index of 
the initial power spectrum of density fluctu ations (n a ) from 
the a nalysis of the 7 year data from WMAP jKomatsu et al.l 
2011). This is the fiducial model and set of priors which will 
also be used in Paper III. For the second analysis, we allow 
for the presence of a non-zero neutrino mass which gives rise 
to a non-zero neutrino density parameter, Q. v . We assume 
uninformative prior information about Q. v (except that it 
is positive-definite) and also continue to use uninformative 
priors on f2 m or erg. For the third analysis, we restore our 
assumption of a negligible amount of neutrino energy den- 
sity, but this time add the parameter wo that describes the 
dark energy equation of state (EoS), such that 



wo = Pa/(c Pa) ■ 



(1) 



Here Pa describes the pressure exerted by dark energy and 
Pa denotes its energy density. 



2.1 The Conditional Luminosity Function 

We use the conditional luminosity function (CLF) model 
described in Paper I to specify our halo occupation distri- 
bution. Readers familiar with our model from Paper I may 
want to proceed directly to Section [3] 



Central CLF 


xLo 


xMi 


71 


72 fc 




9.93 


11.19 


3.5 


0.25 0.156 


Satellite CLF 


6o 


6i 


62 


a s 




-1.08 


1.46 


-0.20 


-1.15 


Nuisance 


'/ 


v> 








1.0 


0.903 







Model parameters used to describe the central and satellite CLFs. 
Bottom row shows the values of the nuisance parameters in our 
model that wc assumed for the Fisher analysis. 



Table 3. Luminosity bins of the SDSS clustering data 



01 M r - 5 log h 


<*) 


«jp accuracy 


AS accuracy 


(1) 


(2) 


(3) 


(4) 


(-19.0, -18.0] 


0.047 


30% 


50% 


(-20.0, -19.0] 


0.071 


30% 


20% 


(-21.0, -20.0] 


0.10 


5% 


15% 


(-21.5, -21.0] 


0.14 


5% 


15% 


(-22.0, -21.5] 


0.17 


5% 


15% 


(-22.5, -22.0] 


0.20 


40% 


15% 



Luminosity bins. Column (1) indicates the magnitude range of 
each luminosity bin (all magnitudes are K+E corrected to z = 
0.1). Column (2) indicates the mean redshift of the lens galaxies 
in each luminosity bin. Column (3) indicates the accuracy of the 
galaxy clustering data and column (4) indicates the accuracy of 
the galaxy-galaxy lensing data. The luminosity function data have 
a fractional accuracy of ~ 4%, and this accuracy degrades to 
about ~ 10 — 15% at the bright and the faint ends respectively. 



The CLF, <$>(L\M) dL, is defined to be the average num- 
ber of galaxies wit h a luminosity L ± dL/2 that reside in a 
halo of mass M (|Yang et all 120031 '). The CLF consists of 
a contribution from central galaxies, $ C (L|M), and from 
satellite galaxies, $ S (L|M). The distribution <& a {L\M) is de- 
scribed by a lognormal distribution with a scatter, <r c , that is 
independent of halo mass, consistent with the findings from 



independent ol halo mass, consi stent with th e finding s trom 
studies of satellite kinemati cs jMore et alj l2009d lb[ |201ll ') 



and galaxy group catalogues l|Yang et alj|2009h . The dep 



dence of the logarithmic mean luminosity, log L c 
mass is given by 



en- 
on halo 



logL c (M) = log 



in 



(M/M1) 71 



[1 + (M/Mi)]' 



C2) 



Four parameters are required to describe this dependence; 
two normalization parameters, Lo and M\ Q and two param- 
eters 71 and 72 that describe the slope of the L C (M) relation 
at the low mass end and the high mass end, respectively. 
The satellite CLF, $ S (L|M) is assumed to be a 



2 The x in front of the Lo and Mi in Table [2] indicates that the 
parameters we use in practice are 10-bascd logarithm of Lq and 



Mi, respectively. 
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Schechter-like function, 



$ s (L|M)dL = $* 



L, 



exp 



AL 
TZ' 



(3) 



Here L* (M) determines the knee of the satellite CLF and is 
assumed to be a factor / s times fainter than L C (M). Moti- 
vated by results from the SDSS group catalog of lYang et al.l 
(|2008l ), we set / s = 0.562 and assume that the faint-end 
slope of the satellite CLF is independent of halo mass. The 
logarithm of the normalization, <3>g is assumed to have a 
quadratic dependence on log M described by three free pa- 
rameters, 6o, bi and 62; 



log $ s * = b + bi (log M - 12) + b 2 (log M - 12) 2 



(4) 



Note that this functional form does not have a physical 
motivation; it merely p rovides an adequa te description of 
the results obtained by lYang et al] (|2008l) from the SDSS 
galaxy group catalog. In addition to specifying the luminos- 
ity dependence of the halo occupation distribution, we also 
need to specify the spatial distribution of galaxies in dark 
matter haloes. Throughout, we assume that central galax- 
ies reside at the center of their haloes and that the satellite 
galaxies follow the matter distribution without any spatial 
bias. T he matter density di stribution is given by the NFW 
profile (|Navarro et al.|[l997h with a concentration that de- 
pends upon mass ac cording to the calibration presented in 
M accio et al.1 |2007t ). As described in Paper I, we allow a 
10 percent uncertainty in the calibration of the normaliza- 
tion of the concentration-mass relation to account for our 
neglect of the baryonic effects on the matter distribution 
of halos and the scatter in the concentration-mass relation. 
We express this uncertainty as a multiplicative parameter 77 
to the fiducial concentration-mass relation and marginalize 
over this nuisance parameter. We use the parameters listed 
in Table [2] for all the three cosmological models. 

In all our models the abundance and clustering of dark 
matter haloes is set by the cosmological parameters, while 
the halo occupation distribution, as parameterized by the 
CLF, describes how this abundance and clustering of haloes 
translates into the abundance and clustering of galaxies as 
well as their cross correlation with matter. In what follows, 
we briefly describe the expressions that can be used to com- 
pute the model predictions for the data described above. 



2.2 Galaxy Luminosity Function 

Given the cosmological parameters and the parameters of 
the CLF, the luminosity function of galaxies, $(!/, z), simply 
follows from multiplying the average number of galaxies in 
a halo of given mass with the number densities of haloes of 
that mass, n(M, z)dM, and by integrating this product over 
all halo masses, 



$(L|M) n(M,z)AM, 



(5) 



[$ C (L|M) + $ S (L|M)] n(M, z) AM . (6) 



It is clear from the above equation, that all of the cosmolog- 
ical information in the galaxy luminosity function is due to 
its dependence on the halo mass function. 

The fraction of galaxies of a particular luminosity 



that are satellites is important to understand the relative 
strength of the luminosity function to constrain the central 
CLF parameters compared to that of the satellite CLF pa- 
rameters. The satellite fraction can be obtained in the CLF 
formalism using 



/sat (L,z) = 



$(L, z) 



$ s (L\M)n(M,z)AM . (7) 



The satellite fraction is also very crucial in determining 
the shape of th e galaxy-galaxy and galaxy-matter cluster- 
ing signals (e gJSeliak et al]|2005l : iMandelbaum et alj|2006l ; 
IZehavi et al.ll201ll) . Typically, the central galaxies dominate 
the galaxy population at all luminosities under considera- 
tion. For the fiducial parameters that we adopt for our analy- 
sis, the satellite fraction is very low at the bright end, but in- 
creases to 30 — 40% at the faint e nd, in good agreement with 
observational constraints (e.g. 



Mandelbaum et al.1 120061 ; 



Ivan den Bosch et al.ll2007l ; iTinker et al.ll2007h 

For the purpose of computing the galaxy-galaxy clus- 
tering and the galaxy-galaxy lensing signal, we will be con- 
cerned with galaxies in a specific luminosity interval [L% , L2] ■ 
The average number density of such galaxies follows from the 
CLF according to 



n s (z) = (N g \M)n(M,z)AM , 



where 



(N S \M) = / $(L|M)dL, 



(8) 



(0) 



is the average number of galaxies with L\ < L < L2 that 
reside in a halo of mass M . 



2.3 Galaxy Clustering 

The clustering of galaxies in the SDSS data is measured us- 
ing a volume limited sample of galaxies and is expressed in 
terms of the galaxy-galaxy correlation function. The corre- 
lation function expresses the excess probability over random 
to find a pair of galaxies separated by a given distance. The 
galaxy-galaxy correlation function consists of two different 
terms based upon the kind of galaxy pairs under consid- 
eration. The correlation function due to pairs of galaxies 
that reside within the same dark matter halo is called the 
one-halo term. On the other hand, the correlation function 
due to pairs of galaxies that reside in separate dark matter 
haloes is called the two-halo term. 

The one-halo correlation can be further subdivided into 
the central-satellite term and satellite-satellite term based 
on the kind of galaxies that constitute the pair. Similarly 
the two-halo correlation can be subdivided into the central- 
central term, the central-satellite term and the satellite- 
satellite term. For computational simplicity, each of these 
correlation function terms are computed in Fourier space. 
The power spectrum and the correlation function form a 
Fourier transform pair. 

The galaxy-galaxy power spectrum, P m (k, z) can be ex- 
pressed as the sum of the following terms, 



P gg (k,z) = 2P c f(k,z) + P s f(k,z) 



+P^(k, z) + 2 P^(k, z) + P^ 



(k, 



(10) 
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As shown in Paper I, these terms can be written in a compact 
form as 

P£(k,z)= f Ux(k,M,z)U y {k,M,z)n(M,z)dM, (11) 



■ (k,z) = J dM 1 H*(k,M 1 ,z)n(M 1 ,z) 

I 



x I dM 2 n y (k,M 2 ,z)n(M 2 ,z)Q(k\M 1 ,M 2 ,z), (12) 

where 'x' and 'y' are either 'c' (for central) or 's' (for satel- 
lite), Q(k\M\, M2, z) describes the power-spectrum of haloes 
of masses Mi and M%, and we have defined 



U c (k, M, z) = H C (M, z) = (iVc [ A f 



and 



Hs(k,M, z) 



{NAM) 
n e (z) 



u s (k\M,z) 



(13) 



(14) 



Here (N C \M) and (N B \M) are the average number of central 
and satellite galaxies in a halo of mass M, which follow from 
Eq. © upon replacing $(L|M) by $ C (L|M) and $ S (L|M), 
respectively. Furthermore, u s (k\M) is the Fourier transform 
of the normalized number density distribution of satellite 
galaxies that reside in a halo of mass M. 

The power spectrum of haloes Q(k\M\, M 2 , 2) is the 
most uncertain component of the halo model, as it needs 
to account for the large scale bias of haloes and the cor- 
responding scale dependence as well as halo exclusion (i.e., 
the fact that haloes are spatially mutually exclusive). We 
do not describe the details of our treatment here but sim- 
ply refer the interested reader to Paper I. We emphasize, 
though, that our analytical model for Q(k\Mi, M2, z) has 
been tested and calibrated using high-resolution numerical 
simulations. Nevertheless, to account for uncertainties, we 
promote one of the calibration parameters, ip, which is used 
to characterize the scale dependence of the halo bias, to be 
a part of our parameter set. Throughout we adopt a 15 per- 
cent uncertainty on tj). 

The corresponding galaxy-galaxy correlation function, 
£ gg (r, z), can be obtained via a Fourier transform of the 
galaxy-galaxy power spectrum, 



1 

2^2 



r{ k,z)^^k 2 dk. (15) 



However, £ gg (j", z) cannot be measured directly from obser- 
vations. This is a direct consequence of our inability to infer 
the line-of-sight separation of galaxies from redshift surveys 
due to the peculiar motions of galaxies. Instead, the correla- 
tion function of galaxies is measured separately as a function 
of the separation of galaxies along the line-of-sight and in 
the plane of the sky to obtain the redshift space correlation 
function Cgg( r pi r Tr,z). This function is then integrated along 
the line-of-sight to obtain the projected correlation function 
w p (r p ,z), 



w p {r p ,z) 



^ s (r p ,r^,z)dr^ . 



(16) 



Here r max is the maximum line-of-sight distance to which 
the redshift space correlation function is integrated in or- 
der to obtain w p (r p ,z). Only for r max = 00, the projected 
correlation function is entirely independent of redshift space 



distortions, and can therefore be expressed in terms of the 
real space correlation function according to 

2 r dr 



w p (r p ,z) = / £ ss (r,z) 



(17) 



However, since real data sets are always limited in extent, 
in practice the projected correlation function w p (r p ,z) is 
always obtained by integrating ^g g (r p ,r 7T ) ou t to some fi- 
nite r max rather than to infinity. For example, IZehavi et al.l 
(2011). whose data we use in Paper III, adopt r max — 
40 /i _1 Mpc or 60 h~ 1 Mpc, depending on the luminosity 
sample used. As discussed in paper I, such values of r max 
are sufficient to get rid of the small scale redshift space dis- 
tortions (commonly called the Finger-of-god effect) but the 
data suffer from residual redshift space distortions on large 
scales that, unless corrected for, can easily re sult in system- 
atic e rrors of 10 percent or larger (s ee also iNorberg et al.l 
120091 : iBaldauf et aT1l2010l : iMorefeoilT ). In Paper I, we have 
shown th at a slightly m odified version of the model pre- 
sented by iKaiserl (| 19871 ) can be used to correct for these 
residual redshift space distortions. Tests using mock galaxy 
redshift surveys show that this method is accurate at the 
few percent level. We refer the interested reader to Paper I 
for further details. 



2.4 Galaxy- Galaxy Lensing 

The clustering of matter around galaxies (i.e., the galaxy- 
matter cross correlation) can be probed using galaxy-galaxy 
lensing; the small distortions of the shapes of background 
galaxies (the sources) due to weak gravitational lensing by 
the matter surrounding foreground galaxies (the lenses). 
Since these shape distortions are extremely weak, back- 
ground galaxies have non-zero ellipticities, and one typically 
can only identify a few background galaxies per lens galaxy, 
a large number of lens galaxies need to be stacked in order 
to detect a signal. The average tangential ellipticity, (e t ), 
around a stack of foreground galaxies is equal to the tan- 
gential shear, jt, which is related to the projected surface 
density around the lens galaxies according to 



AX(R,z) _ E{R,z] - E(R,z) 



{e t )(R,z)= lt (R,z) 

<^crit ^crit 

(18) 

Here, £(i?, z) is the projected matter density averaged 
within a circular aperture of radius R centered around the 
foreground lens galaxy at redshift z and E(_R, z) is the pro- 
jected surface density at a distance R from the lens galaxy. 
The critical surface density £ cl it is a geometric factor that 
includes the angular diameter distances from the observer 
to the source, the observer to the lens and the lens to the 
source. The projected matter density around the foreground 
lens galaxy can be calculated from the galaxy-matter cross 
correlation, £ gm (r, z), using 

f°° 2rdr 

Z(R,z)= p[l + £ gm (r,z)) r—^ . (19) 



R? 



For the purpose of calculating AE, one can safely replace 
[1 + £gm] with £ gm in the above equation. 

The galaxy-matter cross correlation function can be 
modelled using the conditional luminosity function in a man- 
ner that is very similar to modelling the galaxy-galaxy corre- 
lation function: the clustering of matter around galaxies can 
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be divided into one-halo and two-halo terms, each of which 
can be further subdivided into central and satellite terms. 
We calculate each of these terms in Fourier space. The total 
galaxy-matter power spectrum is given by the following sum 

Pgm(k, z) = P cm {k, z) + P sm (k, z) + P cjn (k, z) + P sm (k, z) . 

(20) 

The above terms can be calculated using Eqs. (|11[) - (|12[ ). 
where 'x' is 'm' (for matter) and 'y' is either 'c' (for central) 
or 's' (for satellite). For the matter component, we define 



tions. Carrying out the differentiation of ln£, we obtain 



M 

Um(k, M, z) = —— u h (k\M, z) 
p m (z) 



(21) 



where Uh{k\M, z) is the Fourier transform of the normal- 
ized density distribution of matter within a halo of mass 
M and Pm(z) denotes the average comoving density of the 
Universe at redshift z. The expressions for H c (k, M, z) and 
H a (k, M, z) are given by Eqs. (|13[) and (|14[) . respectively. 
The galaxy-matter correlation function can be obtained by 
a Fourier transform of the power spectrum given by Eq. (|20|) . 
The galaxy-galaxy lensing signal, in turn, can be calculated 
using Eqs. and ([19)) . 

The analytical model sketched above can be used to 
predict the luminosity function of galaxies, their clustering 
strength and the galaxy-galaxy lensing signal around them, 
given a set of cosmological parameters and CLF parame- 
ters. These parameters can be constrained by comparing 
the model prediction to observational data. In what follows, 
we describe the sensitivity of these observations to various 
parameters of our model, and forecast the accuracy with 
which existing data from the SDSS can constrain our model 
parameters, in particular those which describe the cosmol- 
ogy- 



3 FISHER INFORMATION MATRIX 

All the observables in the datasets described above can be 
packed in a single data vector denoted by x. Our model con- 
sists of parameters that describe how galaxies populate dark 
matter haloes, in addition to the cosmological parameters 
that describe the statistics of these dark matter haloes. We 
use the vector to represent the set of our model parame- 
ters. Under the assumption that the probability distribution 
of the data is a Gaussian, the likelihood, £., of the data is 
given by 

In £ = - i In | C| - i (x - M ) T C- 1 (x - M ) . (22) 

Here fi denotes the model predictions at the true parameter 
value and C denotes the covariance matrix of the observ- 
ables. The Fisher information matrix F is given by 



c> 2 ln£ 
dOidQ 



J e 



(23) 



with the derivatives evaluated at and the angular brackets 
denoting an ensemble average over all possible data realiza- 



1 

~2 

1 

2 



d 

do. 



9/i : 

Wi 



— C (x- M ) + (x- M ) C — 



(24) 



C- 1 ((x)- M ) + ((x) 



M) T C" 



dOidOj 



<96>, 89, 



(25) 



(26) 



The terms in the first square bracket are zero because (x) = 
/x(0) and the terms in the second bracket are equal because 
the covariance matrix is symmetric. Therefore, we obtain 

7-1 _ dn T J dn 

no. m, ■ 

The inverse of the Fisher matrix represents the covariance of 
the posterior probability distribution of the model parame- 
ters attainable given the errorbars on the observables in the 
high signal-to-noise ratio (SNR) limit (see lVal lisncri 20081 for 
a more detailed discussion). In particular, the Cramer- Rao 
inequality states that an unbiased estimate of the parameter 
0i given the data can be obtained with an errorbar o~i such 
that 



(27) 



when marginalized over the rest of the parameter set and 
where the equality is satisfied in the high SNR limit. 

For simplicity, we work with dimensionless model pa- 
rameters, Xi, defined by 



= A, I 



(28) 



where 0i denotes the maximum likelihood value for the pa- 
rameter 0i. The Fisher information matrix calculated using 
these new variables is denoted by Fij and is related to Fij 
according to 



p. 



(29) 



With this change of variables, the Fisher information matrix, 
F is also dimensionless and the accuracy with which A; can 
be constrained denotes the fractional accuracy with which 
the parameter can be constrained given the data. In case 
of a diagonal covariance matrix, the dimensionless Fisher 
matrix is given by 



(30) 



where the summation is over all data points. 

As is evident from Eq. (|30[) , the Fisher information ma- 
trix is a (weighted) summation over the derivatives of the 
model predictions, fi, with respect to the model parame- 
ters, 0. The constraining power of any given datapoint on 
a particular parameter depends on the ratio of the corre- 
sponding derivative to the errorbar on this datapoint. The 
larger the absolute value of this ratio, the greater the con- 
straining power. Since we assume that the errorbars are a 
certain fixed percentage of the observables themselves (see 
Table |3]), we are interested in the logarithmic derivatives of 
the observables with respect to the model parameters. These 
derivatives give insight into the power with which each ob- 
servable is able to constrain the model parameter of interest. 
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No priors (n b ,n.,h,Q ln ,a 8 ) WMAP7 priors (n b ,n.,h) WMAP7 priors (D b ,n,,h,0 ln ,(7 8 ) 




Figure 1. Cosmological constraints in the Q m — ag plane that can be obtained by analysing available SDSS measurements of the 
luminosity function (<!>), galaxy-galaxy lensing signal (AS) and the projected clustering measurements (w p ). The 68, 95 and 99 percent 
confidence levels shown by green, yellow and blue contours correspond to the constraints possible when only the <E> and AS data are 
analysed, when only <E> and w p data are used, and when all three $, AS and w p data are used in conjunction, respectively. The three 
panels correspond to different priors assumed on the cosmological parameters as indicated at the top of each panel. Note that, in addition 
to all the CLF parameters and other secondary cosmological parameters, we have also marginalized over the nuisance parameters of our 
model, rj and ?/>. The dotted contours in the middle panel show the 68, 95 and 99 percent confidence levels obtained by WMAP7. 




The logarithmic derivatives of the luminosity function, pro- 
jected galaxy clustering and the galaxy-galaxy lensing signal 
are presented in Appendix A, together with a detailed dis- 
cussion. 



Throughout this paper, we assume that the covariance 
matrix for the luminosity function and the galaxy-galaxy 
lensing signal is diagonal and calculate the Fisher informa- 
tion matrix by using Equation 1301 For the galaxy clustering 
data, we assume that the errorbars are correlated in a man- 
ner which is quantitatively equal to the correlations that ex- 
ist in the measurement s of the projected c lustering of SDSS 
galaxies carried out bv lZehavi et alj (|201lh . We make use of 
Eq. (|29[) to calculate the Fisher information matrix for the 
galaxy clustering data. 



4 PARAMETER CONSTRAINTS AND 
COVARIANCES 

In this section, we forecast the accuracies with which con- 
straints on cosmological parameters can be obtained given 
the accuracy of the current datasets. To obtain these 
bounds, we first calculate the Fisher information matrix, F%j , 
by varying the parameters listed in Tables [T] and [2] We calcu- 
late the Fisher matrix separately for each of the datasets. As 
each of the datasets is independently measured, the Fisher 
information matrix is additive. The inverse of the Fisher 
matrix gives the covariance matrix, C, and the diagonal el- 
ements of this covariance matrix, Cu, represents the accu- 
racy with which the i-th parameter can be constrained after 
marginalizing over the other parameters. In Appendix B, 
we present the procedure we use to obtain 68, 95 and 99 
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Figure 3. Expected cross-correlation coefficients between different model parameters when the luminosity function, galaxy clustering and 
galaxy-galaxy lensing signal are analysed together. Red corresponds to perfectly correlated constraints on the corresponding parameters, 
while blue corresponds to perfectly anti-correlated constraints. The number on the diagonal shows the percentage accuracy with which 
the corresponding parameter can be constrained. 



percent confidence ellipses in a plane corresponding to two 
parameters. In addition to these constraints, we also fore- 
cast the cross-correlation coefficients, Qi j , that are expected 
between the different parameters of our model. We define 
this coefficient according to the following equation, 



\JdiCjj 

The cross-correlation coefficient is a number that ranges 
from [-1,1] and captures the degeneracies inherent in the 
determination of the parameters from the dataset. A posi- 
tive (negative) cross-correlation coefficient, Qij, implies that 
the i-th and j-th parameters are degenerate, such that the 
effect of increasing one parameter can be compensated by 
increasing (decreasing) the other, after allowing the other 
parameters to readjust to the change (which is the essence 
of marginalization). A value close to zero implies that the 
corresponding parameters are only weakly correlated. 

4.1 Model A: Vanilla ACDM 

We first consider the vanilla ACDM model, in which the 
Universe is assumed to have a flat geometry, neutrino mass 
is assumed to be negligible, the initial power spectrum is 
assumed to be a single power-law, and dark energy is mod- 
elled as Einstein's cosmological constant (i.e., has an EoS 



with parameter wo = —1). In order to assess the impact 
of priors on the cosmological constraints, we perform three 
sets of analyses. For the first set, we assume non-informative 
priors for the CLF parameters and all of the cosmological 
parameters. For the second set, we include prior information 
about the secondary cosmological parameters fit,, n s and h 
from the seven year analysis of the cosmic microw ave back- 
ground data from WMAP (jKomatsu et al.ll201ll . hereafter 
WMAP7). For the third set, we additionally include priors 
on the cosmological parameters fi m and as from WMAP7. 
We denote the second set of priors to be the fiducial for 
model A and note that this set will be used in Paper III to 
obtain the constraints from actual data. 

The constraints on the cosmological parameters fi m and 
us resulting from the three different sets of priors are shown 
in Fig. [T] after marginalizing over all CLF parameters, two 
nuisance parameters, r\ and ip, and the secondary cosmolog- 
ical parameters, fib, n s and h. The green contours show the 
68, 95 and 99 percent confidence regions that can be ob- 
tained using a combination of the luminosity function and 
the galaxy-galaxy lensing data. In the absence of prior infor- 
mation on the secondary cosmological parameters, the two 
data sets do not yield particularly tight constraints on fi m 
and as. This is expected in part because the galaxy-galaxy 
lensing signal measured from SDSS lacks information from 
large scales (r p > 2 /i _1 Mpc). As is evident from the mid- 
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Figure 4. Same as Fig. [3] except when only the luminosity 
function and the excess surface density data are analysed. 
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Figure 5. Same as Fig. [3] except when only the luminosity 
function and the projected galaxy clustering data are analysed. 

die panel of Fig. [T] using prior information on the secondary 
cosmological parameters from WMAP7 significantly reduces 
the uncertainties on fi m and erg, although a strong degener- 
acy of the form erg oc fim 3 ' 8 remains. 

The confidence regions highlighted in yellow are a re- 
sult of using a combination of the luminosity function and 
the clustering data. The constraints without the prior infor- 
mation on secondary parameters (left-hand panel) are sig- 
nificantly better than in the previous case. This is primarily 
due to the fact that the clustering data extend out to larger 
radii than the lensing data. As shown in Appendix A, the 
clustering data on large scales (where the 2-halo term dom- 
inates), has good constraining power for fi m and as- Addi- 
tion of the prior information on the secondary cosmological 
parameters from WMAP7 further improves the constraints 



(middle panel of Fig. [TJ . Interestingly, the degeneracies ob- 
tained using the luminosity function and the clustering are 
in an opposite direction to the degeneracies obtained using 
the luminosity function and the galaxy-galaxy lensing sig- 
nal, with as oc Q.^' 2 ■ This immediately suggests that using 
the luminosity function in combination with both clustering 
and lensing will be able to yield even tighter constraints. 

The contours highlighted in blue show the combined 
constraints possible when all three observables are anal- 
ysed together. Even in the absence of informative priors 
on the secondary cosmological parameters, f2 m and as can 
already be constrained to very good accuracy, competitive 
with current constraints in the existing literature (see dis- 
cussion below). In particular, the constraints on f2 m and erg 
from a joint model outperform those obtained by a post- 
combination of constraints from each of the dataset anal- 
ysed separately. A joint analysis is able to effectively con- 
strain the CLF parameters, which is crucial for teasing out 
the tightest possible constraints on the cosmological param- 
eters. This underscores the importance of having each of 
these observables measured from a single uniform survey. 
For comparison, the dotted contours in the middle panel of 
Fig.[T]show the constraints on fl m and erg from the WMAP7 
analysis of the cosmic microwave background anisotropies. 
Note that these constraints are very comparable (in magni- 
tude) to what is achievable using our analysis, albeit with 
roughly orthogonal degeneracies (cf. blue contours in left- 
hand panel). When using the WMAP7 priors on the sec- 
ondary cosmological parameters, the constraints on f2 m and 
erg from the combination of luminosity function, clustering 
and galaxy-galaxy lensing using present-day data from SDSS 
become much tighter than those from WMAP7 alone. In- 
cluding the WMAP7 priors on f2 m and erg results in the 
constraints shown in the right-hand panel of Fig. [1] Because 
the degeneracies on fi m and erg from WMAP7 are roughly 
orthogonal to those resulting from our analysis, the com- 
bination results in extremely tight constraints, even when 
excluding, for example, the lensing or clustering data. 

As indicated earlier, the galaxy-galaxy lensing data that 
we are using to forecast the cosmological constraints only ex- 
tend to projected spatial scales of 2 /i _1 Mpc. In Fig. [2] we 
explore how the constraints on f2 m and erg would improve if 
the galaxy-galaxy lensing was also measured and modelled 
on scales of about 30 h~ 1 Mpc. Such measurements would be 
available in the near future using data from the SDSS (U. 
Seljak, private communication). A comparison of these con- 
straints with those shown in Fig. [2] shows that the largest 
improvement corresponds to the case where we obtain cos- 
mological constraints by combining the luminosity function 
and the galaxy-galaxy lensing data without any priors from 
WMAP7. However, when using the fiducial set of priors on 
the secondary cosmological parameters (or the full set of 
priors from WMAP7), the constraints on £7 m and erg are 
expected to improve only marginally. 

In Fig. [3] we show the cross-correlation coefficients for 
our model parameters expected from a joint analysis of 
the luminosity function, galaxy clustering and galaxy-galaxy 
lensing and using the fiducial prior set. Squares coloured 
in red (blue) indicate the presence of a positive (negative) 
cross-correlation coefficient. Squares devoid of colour indi- 
cate a cross-correlation coefficient close to zero. The number 
on the diagonal squares is the percentage accuracy (rounded 
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to the nearest integer) with which the corresponding param- 
eter can be constrained after marginalizing over all other 
parameters. The first five parameters are the cosmological 
parameters, the next two are the nuisance parameters, fol- 
lowed by the five central CLF parameters and finally the 
4 satellite CLF parameters complete the set. Of the CLF 
parameters, xio, xMi, 72 and a c , which describe the nor- 
malizations, the high mass end slope and the scatter of the 
L C (M) relation are the parameters with constraints that are 
better than ~ 4 percent. The familiar degeneracy between 
xLo and xMi is manifested by a positive cross-correlation 
coefficient (see Appendix A) . In general, the central CLF pa- 



rameters are very tightly coupled with each other and show 
a number of degeneracies. The same is true for the satellite 
CLF parameters. In particular the parameters bo, 61 and 
&2 that describe the normalization, $*(M), are constrained 
very poorly and in a highly correlated fashion. There is also 
a non-negligible cross-talk between the central CLF parame- 
ters, satellite CLF parameters and the nuisance parameters. 

However, this figure clearly shows that the cross- 
correlation matrix nicely separates out into a block diagonal 
form. Most importantly, the block of the cosmological pa- 
rameters is found to have only weak correlations with the 
blocks of CLF parameters and nuisance parameters. This 
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Figure 8. Comparison of the constraints (68 percent confidence) on Q m and erg possible from our joint analysis of the luminosity 
function, galaxy clustering and galaxy-galaxy lensing (shown as a blue shaded region and does not include any prior information on 
the secondary cosmological par ameters) with existing constraints on these parameters from a number of independent m ethods, such as, 
abundance of massive clusters l lVikhlinin et al.ll200gilRozo et alj|2010t iBenson et afll201ll). cosmic shear measurem ent (Hoekstra et al. 
2002l:lLin et al.ll201ll) and halo occupation distribution modelling of galaxy clustering " ijvan den Bosch et al.ll2003l : iTinker et al.ll2005l. 
2012). The WMAP7 constraints are shown using dot-dashed contours. 



implies that the constraints on the cosmological parameters 
are fairly robust to details regarding the modelling of the 
halo occupation statistics. This is in contrast with Fig. [4j 
where we show the cross-correlation matrix expected from 
the combination of the luminosity function and the currently 
existing galaxy-galaxy lensing data. This combination shows 
significant degeneracies between the cosmological parame- 
ters and a number of CLF parameters. These degeneracies 
are suppressed by using the combination of the luminosity 
function and the galaxy clustering data, as can be seen from 
Fig. [5] Analyzing all the three data sets together further 
improves the constraints on fi m and as. 

One of the significant cross-correlation between param- 
eters from the cosmological block and the other parameters, 
when all the three datasets are analysed together, is be- 
tween erg and the nuisance parameter r], which characterizes 
the uncertainty in the calibration of the normalization of 
the concentration-mass relation. This degeneracy is a man- 
ifestation of the well-known fact that the normalization of 
the concentration-mass relation depends on og: dark matter 
haloes in a universe with a larger value of erg collapse ear- 
lier, when the unive rse is denser, and t herefore have higher 
concentrations (e.g.. lMaccio et a"fll2008l ). Hence, an increase 
in erg can be countered by a decrease in 77, which explains 
why the cross correlation coefficient is negative. 

We also expect to obtain excellent constraints on the 
halo occupation distribution of galaxies from our analysis. 
The upper panels of Fig. [6] show the constraints, using the 
fiducial prior set, on the satellite fractions as a function 
of luminosity and the normalization and faint-end slope of 
the satellite CLF as a function of halo mass, respectively, 
for different combinations of data, as indicated in the leg- 
end. The bottom panels show the constraints on the av- 
erage halo mass-luminosity relationship of central galaxies 
and its scatter. It is clear that the combination of luminos- 
ity function and clustering outperforms the combination of 



the luminosity function and the galaxy-galaxy lensing signal 
in constraining the halo occupation distribution parameters, 
and the constraints improve only marginally when analysing 
all three datasets together. The importance of a joint anal- 
ysis, however, is very clear from Figure [7] which shows the 
expected accuracy with which the halo mass function and 
the halo bias function can be constrained. These two quanti- 
ties are the primary source of our cosmological information. 
The constraints on these quantities are fairly broad when 
the galaxy-galaxy lensing data are left out, especially at the 
high mass end. However, they improve considerably with the 
inclusion of the galaxy-galaxy lensing data. 

Finally, in Fig. |HJ we compare the constraints on fi m 
and erg possible from our joint analysis with existing con- 
straints from other independent studies using a variety of 
observation^. To avoid excessive overcrowding in the fig- 
ure, we only show the 68 percent confidence regions from 
each study. As repre sentative of the cosm ic shear studies, 
we show results from lHoekstra et al.l |2002l ) using cyan con- 
tours an d the latest anal ysis of the SDSS Stripe 82 co-added 
data bv lLin et al.l (|201lf ) using black contours. Clearly, the 
current data on cosmic shear measurements is not yet able 
to put interesting constraints on fi m and erg. However, this 
situation is expected to improve drastically in the near fu- 
ture, thanks to a number of deep, large-scale photometric 
surveys planned for this decade. More stringent constraints 
on fim and og have come from studies of cluster abundances. 
As represen tative of these studies , we show the constraints 
obtained bv lVikhlinin et al.l (|2009l ) using X-ray clusters from 
the Chandr a Clus ter Cosmology Project (blue contours), by 
iRozo et al.l (|20ld ) using opt ically selected Max BCG clus- 
ters (green contours), and bv lBenson et alj (|201ll ) using SZ- 

3 The constraints shown in Figure [8] were obtained from 
the respective manusc ripts using the web application Dexter 
jiDemleitner et alJliSplh . 
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selected clusters from the South Pole Telescope (red con- 
tours). 

On this compilation plot, we also show the constraints 
obtained by two studies that used halo occupation mod- 
eling to study the abunda nces and clustering of galaxies, 
van den Bosch et all (|2003r i used the observed abundance 
and luminosity dependence of the correlation length of 2dF- 
GRS galaxies, combined with independent constraints on 
the mass-to-light ratios of galaxy clusters, which resulted 
in the constraints depicted by the br own contours. Using 
a similar method, iTinker et al.l (|2005h modelled the small- 
scale clustering of SDSS galaxies and used the mass-to-light 
ratio on cluster scales to obtain the constraints depicted in 
orange. 



More recently, ITinker et all <|2012h analysed the small 
scale (< 3 /i _1 Mpc) clustering of galaxies from SDSS and 
the mass-to-number ratio of the maxBCG clu ster sample (as 
obtai ned from the weak lensing analysis of ISheldon et al.l 
l2009h to constrain fi m and erg. The 68 percent confidence 
contours from their analysis are shown using magenta con- 
tours. When compared to the forecasted constraints from 
our analysis (shown as a shaded region and including no 
prior information about the secondary cosmological param- 
eters), it is clear that modelling the entire galaxy-galaxy 
lensing signal, rather than just using the mass-to-number 
ratio, yields additional information about the cosmological 
parameters. The fact that the c lustering and lensing sam- 
ples used by ITinker et al.l (|2012h are not well-matched (the 
clustering data is taken from the spectroscopic SDSS, whose 
median redshift is very different from that of the photomet- 
ric sample used to create the maxBCG cluster catalog) also 
introduces additional systematics which result in the weaker 
constraints shown in the figure. 

Note that all these different constraints use different 
priors on the secondary cosmological parameters (or none). 
Therefore their merits cannot be compared directly. How- 
ever, this figure does demonstrate that the joint analysis of 
the abundances, clustering and lensing of galaxies, as advo- 
cated in this paper, is an extremely powerful way of con- 
straining the parameters, f2 m and erg. Such an analysis is 
forecast to yield constraints, even without any additional 
priors on the secondary cosmological parameters, that are 
competitive, if not better than, any of the previous con- 
straints, including the WMAP7 data itself. 



4.2 Model B: Massive Neutrinos 

Next we consider cosmological models with massive neutri- 
nos. We retain the assumption of a flat Universe and the 
standard dark energy EoS, wo = —1. We assume uninfor- 
mative priors on the CLF parameters and the primary cos- 
mological parameters of interest, f2 m , erg and Q. v . The sec- 
ondary cosmological parameters have priors from the 7-year 
WMAP data analysis. The density parameter Q. v is related 
to the mass of the neutrino species, m„ such that 



Q. v h 



94 eV 



(32) 



We perform the Fisher analysis around the cosmological pa- 
rameters for model B displayed in Table [T] The fiducial 
value for the neutrino density parameter is assumed to be 
fi„ = 0.004 which corresponds to ^ m v = 0.184 eV. 



The presence of massive neutrinos in the Universe af- 
fects the matter power spectrum primarily by suppress- 
ing the density fluctuations on scales smaller than the 
neutrino free streamin g scale. We use the framework of 
lEisenstein fc Hul jl999h to obtain the linear power spectrum 
of matter fluctuations in the presence of massive neutrinos. 
iKiakotou et alj (|2008T l have suggested a minor modification 
to this fitting function which leads to a more accurate pre- 
diction of the linear power spectrum when the number of 
neutrino species with degenerate masses N v — 3. We use 
this modification to get the linear power spectrum of matter 
fluctuations. When computing the non-linear matter power 
spectrum, and the halo mass and bias functions from the lin- 
ear theory matter power spectrum, we assume that the cali- 
brations which fit the standard ACDM simulations also gen- 
eralize to the cosmological models that include massive neu- 
trinos. In particular, we assume that the non-linear power 
spectrum is obtained from the linear power spectrum us - 
ing the HALOFIT prescription given bv lSmith et al.l (120031) 
with the modifications suggested recently by Bird et all 
l|2012h based upon numerical simulations including massive 
neutrinos, and that the dependence of the halo mass and bias 
functions on the linear m atter power spectru m is given by 
parameters calibrated bv ITinker et al] (|2008T h In addition, 
we also assume that the density distribution in halos follows 
the NFW profile with the conc entration-mas s relat ion given 
by the calibration obtained bv lMaccio et al] l|2007h . 

The constraints on the parameters fl m , erg and the sum 
of massive neutrinos that can be obtained from a joint analy- 
sis of the abundance, clustering and weak lensing of galaxies 
are shown in Fig. [9] The blue contours represent constraints 
when using the fiducial set of priors on the secondary cos- 
mological parameters. As before, we have marginalized over 
all the CLF parameters, the nuisance parameters and the 
other cosmological parameters. The addition of the param- 
eter, m„ , as expected, leads to more freedom for fi m and 
erg and this slightly weakens the constraints in the Q, m erg 
plane compared to model A. The expected error on the sum 
of neutrinos is not small enough to rule out the massless 
neutrino case (if at least, as assumed here, the true sum of 
neutrino masses is = 0.184 eV). The dotted contours 

show the results from the analysis of the WMAP7 data. It 
is clear that both analyses on their own allow values for the 
sum of neutrino masses as high as 1 eV at 95 percent con- 
fidence. However, the cosmic microwave background data, 
in such a scenario, requires a very low value for erg ~ 0.6. 
Our analysis can rule out such values of erg at very high 
confidence. This shows that combining the constraints from 
our analysis with the cosmic microwave background results 
can certainly provide a significantly better upper limit on 
the sum of neutrino masses. This, in essence, is very similar 
to how the addition of cluster abundances constraints im- 
proves the co nstraints on neutrino masses from the WMAP7 
analysis (e.g.. iBenson et al"1l201ll ). 

The green confidence contours shown in Fig. [9] repre- 
sent the constraints on f2 m , erg and possible from 
our analysis when using all prior information available from 
WMAP7. If the sum of neutrino masses is truly as large as 
we have assumed in our analysis, such a combination will al- 
low the tantalising possibility of the first signs of detection 
of a non-zero sum of neutrino masses at ~ 3cr confidence. 
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Figure 9. The 68, 95 and 99 percent confidence constraints on the parameters ^m„, erg and C m after marginalizing over the CLF 
parameters and other nuisance parameters in our model for two different priors as indicated in the legend. The confidence contours from 
the WMAP7 analysis are shown using dotted lines. 



4.3 Model C: Dark Energy Equation of State 

Finally, we focus on models with a modified EoS for the 
dark energy, i.e., with wo ^ 1. For these models, we restore 
the assumption of massless neutrinos and we maintain the 
assumption of a flat Universe. Modifications to the dark en- 
ergy EoS change the expansion history of the Universe and 
the growth rate of structure formation. The expansion factor 
E(z) is given by 



E{z) = [On(l + zf + Oa(1 + z) 



3(1+™ )i1/2 



(33) 



whereas the dependence of the growth factor on redshift 
can be calculated by solving the second order differential 
equation obtained by combining t he continuity a nd Euler 
equations in the linear regime (see lMo et alll2010h . 



g + 2Hg 



hl m Hl{l + zfg. 



(34) 



Changing the variables from time to f] = lna(t), this equa- 
tion can be written as 



+ 



2 + 



dhiE 



2EH V ) 



(35) 



where a dash indicates a derivative with respect to r\. We 
use the following two boundary conditions to solve for the 
growth factor: (a) g(a = 0.001) = 1 and (b) g'(a = 0.001) = 
1, consistent with the expectation that g oc a in the matter- 
dominated era. We renormalize the scale factor to equal 
unity at redshift z = 0, once the solution has been found. 

These changes to the growth factor that arise from hav- 
ing the dark energy EoS differ from wo = — 1 propagate 
in the redshift evolution of the linear matter power spec- 
trum, and thereby affect the halo mass function, the halo 
bias function and the mass dependence of the concentration 
parameter that describes the density profile of dark mat- 
ter haloes. These changes cause the galaxy abundance, the 
galaxy clustering and the galaxy-galaxy lensing signal to 
deviate from the 'standard' vanilla-ACDM predictions. We 
once again assume that the calibrations required to com- 
pute the non-linear matter power spectrum, the halo mass 
function, the halo bias function and the halo concentrations 
from the linear theory power spectrum are the same as in 
the standard ACDM case. 

The constraints on the dark energy EoS and the cosmo- 
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Figure 10. The 68, 95 and 99 percent confidence constraints on the parameters wo, as, n m from our analysis after marginalizing over 
the CLF parameters and other nuisance parameters in our model for two different priors as indicated in the legend. The confidence 
contours from the WMAP7 analysis are shown using dotted lines. 



logical parameters fi m and as that are achievable from our 
joint analysis, using the WMAP7 priors on the secondary 
cosmological parameters, are shown in Fig. 1101 The dotted 
lines represent the 68, 95 and 99 percent confidence intervals 
from the WMAP7 data, and are shown for comparison. Note 
that our analysis results in degeneracies that are very simi- 
lar to those from the WMAP7 analysis, whereby an increase 
in erg can compensate a decrease in the value of wo- The 
slight differences in the degeneracy directions, however, al- 
low tighter constraints when combining both data sets. The 
resulting confidence intervals are shown using green con- 
tours. 

In Fig. 1111 we compare the constraints on f2 m and wo 
that can be obtained from our joint analysis with those ob- 
tained from other cosmological probes. Supernovae of Type 
la (SNela) and measurements of the baryon acoustic fea- 
ture (BAF) in the clustering of galaxies can be used as 
standard candles and standard rulers, respectively, in order 
to map the expansion history of the Universe as a func- 
tion of redshift. Brown contours indicate the 68, 95 and 99 



percent con fidence intervals on the parameters O m and wo 
obtained bv lSuzuki et~ai] (|2012l ) from the U nion2.1 compi- 
lation of SNela. The constraints obtained bv lPercival et all 
|2010l ). who measured the BAF by combining SDSS and 
2dFGRS data are shown using magenta contours, while the 
constraints from the analysis of the WMAP7 data are shown 
using red contours. Each of these constraints have severe 
degeneracies in the Q m — wo plane, but are very powerful 
when used in conjunction. The constraints forecast for our 
analysis, shown using filled blue contours, indicates that our 
results will be degenerate in yet another direction and can 
therefore significantly add to our knowledge of the cosmo- 
logical parameters f2 m and wq. It is worth mentioning that 
unlike our analysis, the probes of expansion history such as 
SNela are insensitive to the parameter as- The overall am- 
plitude of the clustering measurements used to detect the 
BAF does depend upon as, but cannot be used to constrain 
it, due to its degeneracy with the value of galaxy bias. Our 
analysis, on the other hand, is able to constrain both as and 
galaxy bias simultaneously due to our modelling of small 
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Figure 11. Comparison of the constraints (68, 95 and 99 percent confidence) on Q m and wo possible from our joint analysis of the 
luminosity function, galaxy clustering and galaxy-galaxy lensing (shown as a shaded region and includes WMAP7 priors on secondary 
cosmological parameters) with exi sting constrain t s on these parameters. The magenta contours are from the an alysis of the baryon 
acoustic feature from SDSS data bv lPercival et al l ipoid ), the b rown contours are from the analysis of SNela data bv lSuzuki et al 1 [|2012h 
and the red contours are from the analysis of WMAP7 data bv lKomatsu et ah | toil) . 



scal e clustering and t he galaxy-galaxy lensing signal (see 
also lSeliak et"al1l2005t ). 



5 SUMMARY 

The distribution of galaxies, which is routinely mapped out 
by large scale galaxy surveys, contains a wealth of infor- 
mation about the cosmological parameters. Unfortunately, 
harvesting this information is not straightforward because 
galaxies are biased tracers of the matter distribution. We 
have presented a powerful method that can be used to si- 
multaneously solve for both the cosmological parameters 
as well as the galaxy bias. Our method relies on an accu- 
rate halo occupation distribution modelling of galaxy abun- 
dances, galaxy clustering and galaxy-galaxy lensing, which 
probes scales well into the non-linear regime of structure 
formation. 

In Paper I we developed an analytical framework to 
predict each of these observables, given a cosmology and the 
CLF parameters that describe the halo occupation statistics 
of galaxies. In this paper, we have presented the strength of 
each of these observables (see Appendix A) to constrain our 
model parameters. We have also performed a Fisher matrix 
analysis for a standard ACDM cosmology with flat geometry. 
We primarily focused on the cosmological parameters f2 m 
and as, and treated the other cosmological parameters (fib, 
n s and h) as secondary parameters, for which we adopted 
priors taken from the seven year analysis of the WMAP 
data. 

We have shown that an analysis of the luminosity func- 
tion of galaxies combined with galaxy-galaxy lensing is 
marred by a number of degeneracies between the CLF pa- 
rameters and the cosmological parameters, which effectively 
weakens the cosmological constraints. However, when this 
analysis is complemented with clustering data, a number of 
these degeneracies are broken. For the standard ACDM cos- 
mology, we forecast that the constraints on f2 m and as that 



are achievable from our joint analysis, using data that is al- 
ready available from the SDSS, are of the order of 3 and 2 
percent, respectively. Such constraints will be competitive 
with (and complementary to) the tightest constraints that 
are presently available (see Fig. [5} . We will perform such an 
analysis in Paper III (Cacciato et al, in preparation). 

One of the most important results of this paper is the 
demonstration that the covariance of the posterior distribu- 
tion of parameters from our joint analysis of the luminosity 
function, galaxy clustering and galaxy-galaxy lensing has a 
block diagonal form (see Fig. [3j . In particular, the cosmolog- 
ical parameters form a separate block, which is only weakly 
correlated with the other blocks (2 nuisance parameters, 5 
parameters that describe the CLF of central galaxies, and 4 
parameters that describe the CLF of satellites). This implies 
that the cosmological constraints will be robust to uncertain- 
ties related to details of the halo occupation statistics, and 
that we can simultaneously obtain tight constraints on the 
halo occupation statistics, properly marginalized over the 
uncertainties in the cosmological parameters. 

We have also investigated two extensions of the stan- 
dard ACDM model. In the first extension we included mas- 
sive neutrinos in our analysis. We showed that although the 
constraints on the sum of the neutrino masses, ^2,m v , from 
our analysis are predicted to be of the order of 1 eV (at 95 
percent confidence, assuming realistic errors on the data), 
the degeneracies between rriv and as are oriented differ- 
ently than those from the WMAP7 analysis. Consequently, 
the combination of our results with those from WMAP7 can 
significantly improve the constraints on ^m„. In particu- 
lar, we have shown that such an analysis should be able 
to provide constraints of the order of 0.2 eV at 99 percent 
confidence. 

As a second extension of ACDM, we also allowed for 
a non-standard EoS for the dark energy, characterized by 
wq 7^ —1. We showed that the current data on the abun- 
dances, clustering and lensing of galaxies in SDSS is already 
sufficient to put interesting constraints on wq, especially 
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when combined with the WMAP7 data. In particular, the 
constraints on wo from our analysis are forecast to corre- 
spond to Awq = 0.5 (95 percent confidence, with priors on 
the secondary cosmological parameters). This further im- 
proves to Awo = 0.2 when using the full prior information 
available from WMAP7 (see Fig. [TO]). Since the SDSS data 
considered here only spans a very narrow range in redshift 
(z <, 0.2), there is a huge potential for improvement by mea- 
suring the abundances, clustering and lensing of galaxies at 
higher redshifts. 

Finally, we note that although the main focus of our 
joint analysis is to constrain cosmological parameters, we 
have also shown that we will be able to obtain a detailed, sta- 
tistical description of the galaxy-matter connection, as pa- 
rameterized by the CLF, fully marginalized over uncertain- 
ties in the cosmological parameters. In particular, our anal- 
ysis can provide excellent constraints on the average relation 
between halo mass and central galaxy luminosity, including 
its scatter, on the satellite fraction as function of luminos- 
ity, and on the conditional luminosity function of satellites. 
These constraints characterize, among others, the efficiency 
with which haloes of different masses are able to convert 
their cosmological share of baryons into stars (i.e., galax- 
ies). This will be of great value to inform models of galaxy 
formation and evolution. Another powerful application of a 
tightly constrained CLF is the construction of realistic mock 
galaxy catalogues (see Paper I for an example). The galaxies 
in these mock catalogues will, by construction, have abun- 
dances, clustering and a cross-correlation with matter, that 
are all consistent with current data. 
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APPENDIX A: FISHER INFORMATION 
MATRIX 

The Fisher information matrix, defined as the (negative of 
the) second derivative of the log-likelihood surface, can be 
used to calculate the constraints on the parameters of our 
model, given the observational data. The Fisher information 
matrix can be expressed in terms of the derivatives of the 
observables with respect to the dimensionless parameters Xi 
(see Eqs.[28land [29|) . For our analysis, we have assumed that 
the observational constraints have fixed fractional accuracy. 
Ignoring the covariance between data points for simplicity, 
the dimensionless Fisher information matrix is given by 




where is the fractional accuracy of the k-ih data con- 
straint and [d In /x/t / d\] is the logarithmic derivatives of the 
observable with respect to A 4 . In this appendix, we present 
the logarithmic derivatives of the luminosity function, the 
galaxy-galaxy clustering signal and galaxy-galaxy lensing 
signal with respect to some of our primary model param- 
eters, which enter the Fisher information matrix defined 
above. 

We first focus on the logarithmic derivatives of the lumi- 
nosity function with respect to the central and the satellite 
CLF parameters, shown in different panels of Fig. I All The 
luminosity function of galaxies is always dominated by cen- 
tral galaxies. Therefore it contains more information about 
the central CLF parameters than the satellite CLF param- 
eters, as is evident from a comparison of the magnitudes 
of the logarithmic derivatives in panels (a)-(e) with those 
in panels (f)-(i). If a positive change in a parameter causes 
central galaxies of the same luminosity to live in lower mass 
haloes, which are more abundant, it results in an increase in 
the galaxy luminosity function and it shows up as a positive 
derivative with respect to that parameter. These logarith- 
mic derivatives with respect to the central CLF parameters 
are largest at the bright end (except for the parameter 71 
which controls the low mass end of the halo mass luminos- 
ity relation). On the other hand, the logarithmic derivatives 
with respect to the satellite CLF parameters are largest at 
the faint end because the fractional contribution of satel- 
lite galaxies (satellite fraction) to the luminosity function 
increases as we go to fainter and fainter galaxies. The lumi- 
nosity function gives a large amount of information about 
the parameters Lo and Mi , which determine the pivot points 
of the L c (M) relation. It is interesting to note that that the 
derivatives with respect to xLo and xMi are opposite in 
sign which points to an interesting degeneracy: increasing 
xLo can be compensated by increasing the parameter xMi 
if all other parameters are kept fixed. 

The cosmological parameters affect the luminosity func- 
tion by changing the halo mass function. As can be seen from 
panels (a)-(c) of Fig |A2l most of the information about the 
cosmological parameters comes from the bright end of the 
galaxy luminosity function. This is expected because the 
change in the halo mass function due to the cosmological 
parameters is largest at the most massive end. The deriva- 
tives with respect to fi m and as have the same sign and 
are only slightly different in shape. Hence a positive change 
in f2 m can be countered with a negative change in <rg and 



vice versa. The weakest constraints are on the parameter wo- 
The parameter wq mainly affects the luminosity function by 
changing the growth of structure and thus the power spec- 
trum. We have assumed that the luminosity function is only 
measured at a single redshift z — 0.1 from SDSS data, thus 
the leverage in redshift to propagate changes to the growth 
of structure are not very large. This results in very weak 
sensitivity to wo- 

The logarithmic derivatives of the projected galaxy clus- 
tering with respect to the CLF parameters are shown in 
Fig. IA3I On large scales, the clustering of galaxies mainly 
arises due to the two halo term which is sensitive to the clus- 
tering of haloes in which these galaxies reside. Since more 
massive haloes are more strongly clustered, a change in a pa- 
rameter which causes galaxies to reside on average in more 
massive haloes will result in an increase in the clustering 
of galaxies on large scales. On small scales, the one halo 
terms (central-satellite and the satellite-satellite) dominate 
the clustering of galaxies. The satellite distribution in more 
massive haloes is less concentrated but the number of satel- 
lites in these haloes is larger and there are competing effects 
that can result in changes in the clustering of galaxies in di- 
rections opposite to the trend on large scales (see e.g. panel 
[d]). The trends in the logarithmic derivatives, as shown in 
Fig. IA31 on average show that the small scale clustering of 
galaxies has the most information about the halo occupa- 
tion parameters. The galaxies in the brightest bins show a 
large derivative with respect to the central CLF parame- 
ters (barring 71 as expected). The satellite CLF parameters 
that determine the normalization of the satellite CLF are 
almost equally sensitive to all luminosity bins. The logarith- 
mic derivatives are largest for the parameters xLo and xMi. 
Although the derivatives are opposite in sign, the shapes of 
these derivatives as a function of scale are quite different 
from each other, thus providing an avenue to constrain the 
parameters separately. 

Comparisons of the logarithmic derivatives of the galaxy 
clustering signal with respect to the cosmological parame- 
ters, shown in Fig. IA4I to the derivatives shown in Fig. IA21 
give an interesting insight. Unlike the logarithmic deriva- 
tives of the luminosity function, the derivatives with respect 
to H m and as have opposite signs, which suggests that a 
combination of the two observables adds a lot of informa- 
tion about these parameters. The derivatives of the galaxy 
clustering signal with respect to the parameter fi m increase 
on large scales. The parameter erg controls the concentration 
of dark matter haloes and therefore the distribution of the 
satellite galaxies. Therefore it affects both the small and the 
large scales. Galaxy clustering on large scales is also sensi- 
tive to the parameter wo, because it affects the growth of 
structure, while on small scales the sensitivity is expected 
due to the changes to the concentration-mass relation. 

The information contained in the galaxy-galaxy lensing 
signal, as shown by the logarithmic derivatives in Fig. IA5I 
and Fig. IA61 is dominated by the one halo terms as the sig- 
nal has only been measured out to a couple of megaparsecs. 
The derivatives with respect to the CLF parameters show 
quite different dependencies on scale compared to the clus- 
tering data, with a large amout of information coming from 
the scales of the order of a megaparsec. Most of the trends 
are fairly similar to the trends in the logarithmic derivatives 
of the galaxy-galaxy clustering data and depend upon how 
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Figure Al. The logarithmic derivatives of the luminosity function with respect to the dimensionless CLF parameter A. The CLF 
parameter corresponding to each logarithmic derivative is indicated at the top of each panel. 
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Figure A2. The logarithmic derivatives of the luminosity function with respect to the dimensionless cosmological parameters A. The 
cosmological parameter corresponding to each logarithmic derivative is indicated at the top of each panel. 
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Figure A3. The logarithmic derivatives of the galaxy clustering signal with respect to the dimensionless CLF parameter A. Different 
line types correspond to the six different luminosity bins. The faintest bin is shown using solid line, while the brightest bin is shown 
using a dot-long-dashcd line. The CLF parameter corresponding to each logarithmic derivative is indicated at the top of each panel. 
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Figure A4. The logarithmic derivatives of the galaxy clustering signal with respect to the dimensionless cosmological parameter A. 
The cosmological parameter corresponding to each logarithmic derivative is indicated at the top of each panel. The different line types 
correspond to the same bins as that in Fig [A3] 
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Figure A5. The logarithmic derivatives of the galaxy-galaxy lensing signal with respect to the dimensionless CLF parameter A. The 
CLF parameter corresponding to each logarithmic derivative is indicated at the top of each panel. The different line types correspond to 
the same bins as that in Fig |A3l 
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Figure A6. The logarithmic derivatives of the galaxy-galaxy lensing signal with respect to the dimensionless cosmological parameter A. 
The cosmological parameter corresponding to each logarithmic derivative is indicated at the top of each panel. The different line types 
correspond to the same bins as that in Fig [A3] 
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the CLF parameters change the halo occupation distribution 
of galaxies. The logarithmic derivatives with respect to the 
cosmological parameters show that measuring the galaxy- 
galaxy signal to larger scales can certainly add more infor- 
mation. We have quantified how the constraints on the cos- 
mological parameters improve by using the galaxy-galaxy 
lensing measured to ~ 30 /i _1 Mpc signal in Section [4] 



to the eigenvectors, and with axes lengths equal to 1.5096, 
2.4477 and 3.0349 times the square root of the correspond- 
ing eigenvalues found in step (iii), respectively. Note that, 
this procedure can be trivially generalized to obtain confi- 
dence ellipsoids in a subspace spanned by any number of 
parameters. 



APPENDIX B: CONFIDENCE INTERVALS 
FROM THE FISHER MATRIX 

The inverse of the Fisher information matrix gives the co- 
variance with which the parameters of the model can be con- 
strained given the data. Here, we briefly present the proce- 
dure to obtain the confidence intervals in a two dimensional 
parameter space given this covariance matrix. Let P(x,y) 
denote the posterior probability distribution for two param- 
eters of interest x and y after marginalizing over the rest 
of the parameter set. For simplicity, let us assume that the 
constraints on these two parameters are uncorrelated, which 
implies that the subset of the full covariance matrix corre- 
sponding to these two parameters is a diagonal matrix given 

• .-Ct i). » 

In case the covariance matrix is not diagonal, one can al- 
ways diagonalise it by rotating the co-ordinate system and 
aligning the axes with the eigenvectors of this covariance ma- 
trix. Let us assume that the posterior distribution P(x, y) is 
Gaussian and without loss of generality, let us assume that 
it is centred at (0,0), such that 



P ( x > V) = r. u ~ i ex P 
27r \a xx a yy \ 



(B2) 



Changing variables to the dimensionless quantities £' and 8 
defined such that 

x = a xx £ cos(#); y = a yy ^' sin(0) , (B3) 

gives 

P(x,y)dxdy = P(Z',0)dt'd6, (B4) 



2tt 



exp . (B5) 



The probability that the true model lies within a region 
enclosed by the iso-probability contour £ is then given by 



Pt= I [ P(£',0)d£'d0 = 1 -exp 



-e 1 



(B6) 

The iso-probability contours £ = 1.5096, 2.4477 and 3.0349 
correspond to P% = 0.68, 0.95 and 0.99, respectively. These 
values of £ correspond to 68, 95 and 99 percent confidence 
ellipses in the x — y plane with principal axes of length equal 
to £a xx and £<r y y- 

To summarize, the step-by-step procedure we use to ob- 
tain the confidence ellipses is as follows: (i) invert the entire 
"-par x n par Fisher matrix to obtain the covariance matrix, 
(ii) choose the 2x2 subset of the covariance matrix cor- 
responding to the two parameters of interest, (iii) find the 
eigenvalues and eigenvectors of this submatrix, (iv) draw 68, 
95 and 99 percent confidence ellipses with the axes aligned 
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